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ABSTRACT 

We discuss pulsations of the rapidly oscillating Ap (roAp) star HD 24712 (HR 1217) 
based on nonadiabatic analyses taking into account the effect of dipolc magnetic fields. 
We have found that all the pulsation modes appropriate for HD 24712 are damped; 
i.e., the kappa- mechanism excitation in the hydrogen ionization layers is not strong 
enough to excite high-order p-modes with periods consistent with observed ones, all 
of which are found to be above the acoustic cut-off frequencies of our models. 

The main (2.721 mHz) and the highest (2.806 mHz) frequencies are matched 
with modified / — 2 and / — 3 modes, respectively. The large frequency separation 
(« 68 /iHz) is reproduced by models which lay within the error box of HD 24712 on 
the HR diagram. The nearly equally spaced frequencies of HD 24712 indicate the small 
frequency separation to be as small as « 0.5 /iHz. However, the small separation 
derived from theoretical / = 1 and 2 modes are found to be larger than ~ 3 /zHz. 
The problem of equal spacings could be resolved by assuming that the spacings cor- 
respond to pairs of / = 2 and I = modes; this is possible because magnetic fields 
significantly modify the frequencies of I = modes. The amplitude distribution on 
the stellar surface is strongly affected by the magnetic field resulting in the predom- 
inant concentration at the polar regions. The modified amplitude distribution of a 
quasi-quadrapole mode predicts a rotational amplitude modulation consistent with 
the observed one. 

Amplitudes and phases of radial-velocity variations for various spectral lines are 
converted to relations of amplitude/phase versus optical depth in the atmosphere. 
Oscillation phase delays gradually outward in the outermost layers indicating the 
presence of waves propagating outward. The phase changes steeply around logr ^ 
—3.5, which supports a T — r relation having a small temperature inversion there. 
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1 INTRODUCTION 

HD 24712 (HR 1217, DO Eri) is a prototype of the rapidly 
oscillating Ap (roAp) stars that consist of about 40 mem- 
bers. The oscillations, with periods ranging from ~ 6 to 
^ 20 min, are high radial-order p-modes affected by strong 
magnetic fields of 1 kG to 25 kG. 

The 6. 15 min u te lig h t variation of HD 2471 2 was dis- 
covered by iKurtj (|l98ll ). iMatthewT et all (|l988l ) first dis- 
covered radial velocity (RV) variations with an amplitude of 
400±50ms"^ A WET (Whole Earth Telescope) campaign 
found ei ght frequenc i es wi th rotational side-lobes as pre- 
sented in iKurtz et al.ll|2005h : the paper also gives a thorough 
review of precedi ng research on this star. F rom accurate 
RV measurements lMkrtichian fc Hatzed (|2005l ) obtained two 
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additional frequencies. iRvabchikova et al] (|2007| ) (Paper I) 
studied oscillation amplitudes and phases corresponding to 
the two highest amplitude modes for «600 unblended spec- 
tral lines of different elements/ions, and found a diversity of 
these pulsational characteristics in two thirds of them. 

The magnetic field of HD 24712 has been studied b y 
many authors as re viewed in iRyabchikova et al] l|l997h . 
iBagnulo et all l| 19951 ) derived a polar magnetic strength of 
Bp = 3.9 kG, a rotational inclination (angle between line-of- 
sight and rotation axis) of i = 137° (or 43° ) and a magnetic 
obliquity (the angle bet ween the rotation a nd magnetic axes) 
of /? = 150° (or 30°). Rva bchikova et al.l lfl997. 1 found that 
a slightly higher strength of Bp = 4.4 kG is more consistent 
with their polarimetric observations. Recen tly, by inverting 
rotationally modulated polarized spectra, iLiiftinger et al.l 
l|2008l ) found that HD 24712 has a nearly dipole magnetic 
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field with the strength varying between 2.2 kG and 4.4 kG, 
depending on the rotation phase. 

One of the reasons why oscillations of roAp stars are im- 
portant lies in the fact that many (often regularly spaced) os- 
cillation frequencies are excited simultaneously giving great 
potential for asteroseismic studies. One difficulty in apply- 
ing asteroseismology to roAp stars is that strong magnetic 
fields affect the high-order p-modes in complex ways. Al- 
though the frequency of an oscillation mode generally in- 
creases with increasing magnetic field strength, it occasion- 
ally jumps down by ~ 10 to 30 /xHz and then starts increas - 
ing again (|Cunha fc Goughl [2OO0I : [S aio fc Gautschvll2004l ). 
At present, there are three independent methods to calcu- 
lat e magnetic effects on h igh-order p-modes: the method 
of ICunha fc Goughl l|2000l ) is based on a variational prin- 
ciple, the other two use truncated expansions with spheri- 
cal harmo nics to present the angu l ar de p endences of eigen- 
functions llDziembowski fc Goodd 1 19961 : iBigot et al.1 12OO0I : 
ISaio fc Gautschvll2004l ). The results from those thre e meth- 
ods ro ughly a gree with e ach other as discussed by ICunhal 
(|200d ) and bv lSaid (|2008l '). 

A furthe r complexity arises if th e rotat ion effect is taken 
into account. iBigot fc Dziembowskil l|2002l ') investigated the 
oscillation properties taking into account both effects of ro- 
tation and magnetic field. In our present investigation, how- 
ever, we disregard the effect of rotation on the frequency 
and eigenfunctions for simplicity, hoping that the effect is 
small because the pulsation periods, ~ 6 min, are very much 
shorter than the « 12-d rotation period of HD 24712. 

Previo usly, t heore tical pulsation models based on the 
method of ISai 

d (|2005l ) were compared with the observed 
frequencies for 7 Equ l|Gruberbauer et aL 20081). 10 Ag l 
jHuber et al.ll2008l ) and HD 101065 (|Mkrtichian et al.ll2008f) . 

Those attempts were more-or-less successful, although the 
required strength of the magnetic fields in models tended 
to be larger than that measured by the Zeeman splittings of 
spectral lines. In this paper we provide the same modelling of 
HD 24712 and for the first time we compare not only oscilla- 
tion frequencies but also the amplitude and phase variations 
in the atmosphere of HD 24712. 



2 MODELS 



Table 1. Unperturbed model parameters 



Name 


M/Mq 


X 


Z 


r(r) 


He-dep. 


conv 


AD160 


1.60 


0.70 


0.02 


Ap 


Y 


N 


AD165 


1.65 


0.70 


0.02 


Ap 


Y 


N 


AD170 


1.70 


0.70 


0.02 


Ap 


Y 


N 


AH165 


1.65 


0.70 


0.02 


Ap 


N 


N 


AH165C 


1.65 


0.70 


0.02 


Ap 


N 


Y 


AD170Z25 


1.70 


0.695 


0.025 


Ap 


Y 


N 


AH170Z25 


1.70 


0.695 


0.025 


Ap 


N 


N 


AD150Z1 


1.50 


0.71 


0.01 


Ap 


Y 


N 


AH150Z1 


1.50 


0.71 


0.01 


Ap 


N 


N 


SD160 


1.60 


0.70 


0.02 


ss 


Y 


N 


SD165 


1.65 


0.70 


0.02 


ss 


Y 


N 


SD170 


1.70 


0.70 


0.02 


ss 


Y 


N 


SH165 


1.65 


0.70 


0.02 


ss 


N 


N 



0,7 




-4 -2 

log T 



Figure 1. Two 

standard relation 
IShibahashi fc Sai o' 



types of T 
dash ed line) 

1985'), while the Ap star relation (solid line) 
is adopted from iShulvak ct al.. (.2009.) 



relations employed. The 
IS obtained from eq. (10) of 



2.1 Unperturbed models 

For unperturbed models, we adopted spherically symmet- 
ric e volutionary models calc ulated with OPAL opacity ta- 
bles (llglesias fc Rogers|[l99^ ). Table [1] lists various assump- 
tions adopted for each series of evolutionary models. Except 
for AH165C, envelope convection is suppressed, assuming a 
strong magnetic field to stabilize convection. For the enve- 
lope convection included in AH165C we used a local mixing- 
length theory with a mixing-length of 1.5 times the pressure 
scale height. 

We considered two types of T — r relations in the opti- 
cally thin outer layers (Fig.[TJ. One of them is a standard re- 
lation (denoted as " SS" in Table[T} obtained from eq. (10) of 
IShibahashi fc Said l| 19851 ). The other one (denoted a s "Ap" 
in Table [TJ is a relation obtained by IShulvak" et al.1 (|2009l ) 
for the self-consistent atmospheric model of HD 24712 that 
takes into account element stratification including rare-earth 



elements. The latter relation has a small temperature inver- 
sion at logrsooo ~ — 3.5 where abundances of praseodymium 
(Pr) and neodymium (Nd) increase steeply outward. (Al- 
though the optical depth in Shulyak et al.'s relation refers 
to rsooo, optical depth at a wavelength of 5000 A, we as- 
sume in our present paper that rsooo is not very different 
from Rosseland-mean optical depth r.) 

We have computed models with and without helium 
depletion in the outermost layers (the 6th column of Ta- 
ble [!}. In the helium depleted models, the helium abun- 
dan ce is calculated as Y = 0.01 + {Yi — 0.01) x {x2 + xs) 
fcf. ' Balmforth et al.ll200ll ). where Yi is the (initial) helium 
abundance in the interior, X2 and X3 are fractions of He II 
and He HI, respectively. 

Fig-Elshows the evolutionary tracks of AD160, AD165, 
AD170, AD170Z25, and AD150Z1 with the position of 
HD 24712 with error bars. The type of T — r relation, on/off 
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Figure 2. Evolutionary tracks and tlie position of HD 24712 with 
error bars. Solid lines show the evolutionary tracks for a stan- 
dard composition of {X,Z) = (0.7,0.02). Open squares indicate 
the loci of models giving best fits to the observed frequencies for 
various cases listed in Table [T] At those positions, models have 
large frequency separations similar to 68 fiHz, which is determined 
mainly by stellar mass and radius. 



of the He-depletion, and on/ofT of the envelope convection 
hardly change the evolutionary tracks on the HR diagram. 
The stellar masses for the metal rich {Z = 0.025) and metal 
poor {Z = 0.01) cases are chosen such that the evolutionary 
tracks pass close to the position of HD 24712. 

We have adopted the effective temperature Tag — 
7250 ± 150 K (or l og Teg = 3.860 ± 0.009 ) obtained by 
iRvabchikova et all (|l997l ) and confirmed by IShulvak et al.l 
I 2009h . The range of effective temper ature i s cons istent 
with Teft = 7330 ± 140 K derived by [Wadd ll 19971 ) and 
Tefi = 7350 K dete r mined bv lLiiftinger et al. I ([2008'). 

IShulvak et all (|2009l ) estimated the radius of HD 24712 
to be 7? = 1.772±O.O43R0 by combining the spectral energy 
distribution and the Hipparcos parallax vr = 20.32±0.39 mas 
Jvan Leeuwenll2007[ ). These estimates for the effective tem- 
perature and radius yield logL/L© = 0.891 ± 0.041 for the 
luminosity of HD 24712, which was used to place HD 24712 
on the HR diagram (Fig. [2]). We note that the adopted lu- 
minosity and temperature of HD 24712 are consistent with 
the values used in previous analysis by Cunha ct al. (2003) 
(logL/L0 = 0.892±0.041 and logT^ff = 3.869i!J;;iJ?|), where 



the l uminosity (based on th e Hipparcos parallax) was taken 
from lMatthews et al.l l|l999l ). 



2.2 Pulsation Models 

We have obtained nonadiabatic frequencies and eigenfunc- 
tions for axisymmetri c hig h -order p-modes based on the 
method described in I Said (|2005l ) except for the outer- 
boundary condition and the perturbation of radiation. Since 



all of the observed oscillation frequencies in HD 24712 are 
above the acoustic cut-off frequencies of the models consid- 
ered here, we used a running wave condition for the mechan- 
ical outer boundary condition at log r ~ —6. 

Applying a standard method (see e.g.. lUnno et al.l 1989f) 
to the mechanical equations (A3) and (A4) of ISai o n|2005h . 
we obtain a running-wave condition as 



V 

2Xp 



■i^ALu'^Xp/V 



(1) 



where Xp = {din P/d\n p)t with P and p being pressure 
and matter density, respe ctively, and the other symbols are 
the same as in'Saio (2005). In deriving the equation we have 
assumed the oscillations are isothermal at the outer bound- 
ary. 

In calculatin g the perturbat i on of radiation flux, we 
have adopted the lUnno fc Spiegell (|l966l ) theory for the Ed- 
dington approximation, where the radiative flux is given as 



3kp^ ( 



acT* 



ldS_\ 
n dt J ' 



(2) 



where S, k, a, and c are entropy per unit mass, opacity per 
unit mass, the radiation constant, and the speed of light, 
respe ctively. We used a linearized form of the above equa- 
tion (ISaio fc Coxlll980l ). Including the dS/dt term with the 
running wave boundary condition is found to reduce the os- 
cillation phase variation in the outermost layers making it 
consistent with observations. 

The angular dependences of eigenfunctions are ex- 
panded using axisymmetric spherical harmonics Y^™^" (i.e., 
Legendre polynomials) with i = 2j — 2 (even modes) or 
£ = 2j — 1 (odd modes), where j — l,2,...jt. An even 
(odd) mode is symmetric (anti-symmetric) with respect to 
the magnetic equator. The expansion is truncated as jt = 12 
in most cases; but sometimes jt = 14 is adopted to have bet- 
ter accuracy in frequencies. To label the type of angular de- 
pendence of a mode we use Im which is equal to the £ value of 
the component having the maximum kinetic energy among 
the expansion components. We note that the angular dis- 
tribution of amplitude varies depending on the strength of 
magnetic fleld even for a flxed Im ■ To represent the strength 
of a dipole magnetic field we use the polar strength Bp. 



3 THEORETICAL RESULTS AND 

COMPARISONS TO OBSERVATIONS 

3.1 Oscillation frequencies 

We calculated nonadiabatic oscillation frequencies between 
~ 2.5 mHz and ~ 2.9 mHz for ^ l,n ^ 3 including the ef- 
fect of dipole magnetic fields in a range of 2 Bp/ (kG) ^ 7. 
For the series listed in Table [U we performed pulsation anal- 
yses for models having large frequency spacings comparable 
with that of HD 24712; i.e., 68 /iHz. 

All the pulsation modes having frequencies compa- 
rable with the observed ones are found to be damped; 
i.e., no excited modes are found. The excitation mecha- 
nism for the oscillations of roAp stars is generally thought 
to be the K-mechanism in the hyd rogen ionization zone 
(|Balmforth et all l200ll : ICunhal |2002|). For cool roAp star s 
like HD 24712 and HD 101065 (|Mkrtichian et al.l |2008| ). 
however, the K-mechanism seems not strong enough. The 
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Figure 3. Observed oscillation frequencies and amplitudes of 
HD 24712 (left panel) and theoretical frequencies as a function 
of Bp (the strength of magnetic field at poles) for a model in the 
series AD165 (right panel). In the left panel, solid lines show fre- 
quencies and the photometric amplitudes obtained by the WET 
campaign (fl-f8; iKurtz et al1l2005, '). while dotted lines represent 
additional two frequenc ies (fml and fm2) obt a ined by the radial 
velocity observation by iMkrtichian &: HatzesI l(200^ , the ampli- 
tudes of which are set arbitrary equal to the amplitude of the 
WET lowest frequency. The label for each frequency is indicated 
along the rightmost vertical axis. The middle panel shows theoret- 
ical frequencies obtained without including the magnetic effects. 
Horizontal dotted lines in the right panel indicate the observed 
frequencies. 



temperature inversion (Fig. [1} is t oo small to help excit e 
high-order p-modes as discussed in iGautschv et al.l l)l998h . 
It seems that we need a new excitation mechanism for cool 
roAp stars. 

Fig. [3] shows an example of p-mode frequencies in a 
model as a function of Bp (the magnetic-field strength 
at poles). The obs erved frequencies by a WET campaign 
IIKurtz et al. I2OO5I ) and the radial velocity observation by 
IMkrtichian et alTlj2008l ) are also shown. Different symbols 
correspond to different values of Im- Generally, the frequency 
of an oscillation mode increases with Bp, but occasionally it 
jumps down by 30 pHz and the n starts increasing again . 
This property was discovered by ICunha fc Goughl (I2OO0I ) 
and confirmed by I Saio fc Gautschvl (|2004l ). The expansion 
method adopted in this investigation fails around the fre- 
quency jump, where the distribution of kinetic energy among 
the components of expansion is broad. 

Frequencies of different radial order for a given Im vary 
approximately parallel as a function of Bp so that the large 
frequency separation Ai/ hardly changes with Bp, where 
Au = u{lm',n) — u{lm',n — 1) with n being the radial or- 
der of a mode. The top panel of Fig. 3] shows Av calculated 
for Im = 2 modes in the model shown in Fig. |3] The large 
separation increases very slowly as a function of Bp. This 



S = v(l = \,n) - 0.5*[L'(J=2,n-l) + y(f=2,n)] 




Hp(KG) 



Figure 4. Large separation (top panel) and small separation 
(bottom panel) of frequencies as functions of Bp for the model 
shown in Fig. [3] Dotted lines indicate observed large and small 
separations of HD 24712. 



corresponds to an increase in the phase velocity of magneto- 
acoustic wave, yj ci + v\, with Bp, where Ca and va are the 
adiabatic sound speed and the Alfven speed, respectively. 

The bottom panel of Fig. |4] shows the small separation 
defined as 

5 = u{lm = 1; n) - Q.^[v{lm = 2; n - 1) + u{lm = 2; n)] (3) 

calculated using the frequencies shown in Fig. [S] The small 
separation of the model is about 5 /iHz at Bp — 0. In the 
presence of a magnetic field, 5 varies from ~ 3 fiRz to ~ 
6 //Hz. Generally, 5 increases with Bp, but it jumps down at 
the frequency jumps discussed above. The gradual growth 
of 5 comes from the fact that the frequencies of = 1 
modes increase slightly more steeply with Bp than those of 
Im = 2 modes. Just after the jump, S reaches a minimum 
which is considerably smaller than in the case of Bp — 0, but 
still significantly larger than the observed value 0.5 /iHz, 
obtained by assuming that f3 is a = 1 mode and f2 and 
f4 are = 2 modes. Such a small value of S cannot be 
reproduced by any model examined in the present paper. 
The problem of the small spacing of HD 24712 is not unique. 
iBruntt et all (|2009l ) found a small spacing that is essentially 
zero in another roAp star a Cir - very much smaller than 
that of HD 24712. It is important to solve this problem of the 
small spacing; it could be related to a fundamental property 
of roAp stars. 

A relatively large theoretical 5 means that the frequency 
of an Im = 1 mode is slightly larger than the mean of the 
adjacent two Im = 2 modes. It is interesting to note that 
Im ~ modes cross Im = 1 modes at Bp « 5 kG (Fig. |3]). In 
other words, the frequency of each — mode is slightly 
smaller than the frequency of the adjacent Im = 1 mode 
at Bp « 5 kG. Therefore, the problem of equal spacings of 
HD 24712 can be solved if we assign Im = modes (rather 
Im = 1 modes) to fl, f3, and f5, and consider that Bp hap- 
pens to have a value near the crossings between Im = and 
Im ~ i modes. The solution has some weaknesses, however: 




0.26 0.27 



log R 

Figure 5. The distribution of the mean deviation of theoretical 
frequencies from the observed ones on the plane of Bp and stellar 
radius along the evolution sequence AD165. Darker parts have 
smaller mean deviations. Contours drawn for mean deviations of 
1.5 /iHz, 2 ^Hz, and 3 /jHz. The white areas have mean deviations 
greater than 10 /iHz. Although the mean deviation in this diagram 
was calculated without including Im = modes, the distribution 
of the mean deviation with 1^ = modes is similar. 



1) For the Im = (2, 0, 2) combination to become equally 
spaced, Bp should have a particular value at which the fre- 
quency of an Im ~ mode is just equal to the mean of the 
two adjacent Im = 2 modes; it looks unlikely for such a co- 
incidence to be realized in a Cir, too; 2) As we will discuss 
below, the amplitude of an = mode modulates with 
rotation phase differently from the observed total hght vari- 
ations. Nevertheless, we will consider both possibilities for 
the mode identifications for fl, f3, and f5. 

The observed frequencies of HD 24712 were compared 
to models having Av ~ 68 fiRz for each evolution sequence 
listed in Table [T] The quality of the fits was estimated by 
the mean deviation of model frequencies (for Im ^ 3) from 
the nine observed ones. The WET f7 frequency was not in- 
cluded in calculating the mean deviation because it differs 
from f6 only by 2.4 /xHz and no models can fit both fre- 
quencies. In calculating mean deviations, we adopted equal 
weights for all the frequencies because observational errors 
are much smaller than the theoretical uncertainties which 
can be estimated as ~ 1 /iHz from the scatter seen in Fig. 

Table [2] lists parameters and the mean deviation (MD) 
of the best-fit model for each model series, where MDs in 
parentheses refer to the values obtained including Im = 
modes. Generally, including Im = models yields much bet- 
ter agreement. The positions of these best models on the 
HR diagram are shown by square symbols in Fig. (2] The 
position of HD 24712 on the HR diagram is consistent with 
1.65 M© model with normal composition and with 1.70 M0 
model with a heavy-element abundance oi Z = 0.025. The 
quality of the frequency fit is practically independent of the 
type of T-T relations, helium depletion, or the efficiency of 
convection. The mean deviations of the best-fit models tend 
to be smaller when the position of the models on the HR di- 
agram is close to the position spectroscopically determined 
for HD 24712.^ 

Ahhough ICunha et all (|2003l ) concluded that a metal- 
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Table 2. Best fit models 



Name 


logR 


lo 




logL 


Bp 


MD(^tHz) 






Q 
O 


8471 
1 ± 




o.o 


1.54 


(1.09) 


AD165 


0.2690 


3 


8585 


0.9247 


4.9 


1.45 


(1.04) 


AD170 


0.2724 


3 


8696 


0.9759 


4.6 


1.60 


(1.05) 


AH165 


0.2702 


3 


8581 


0.9255 


5.5 


1.45 


(1.03) 


AH165C 


0.2723 


3 


8574 


0.9267 


6.8 


1.40 


(1.06) 


AD170Z25 


0.2735 


3 


8549 


0.9192 


5.1 


1.48 


(0.99) 


AH170Z25 


0.2668 


3 


8572 


0.9150 


5.6 


1.45 


(1.43) 


AD150Z1 


0.2548 


3 


8667 


0.9291 


4.6 


1.71 


(1.22) 


AH150Z1 


0.2561 


3 


8663 


0.9299 


5.3 


1.64 


(1.12) 


SD160 


0.2665 


3 


8467 


0.8726 


5.8 


1.49 


(1.06) 


SD165 


0.2700 


3 


8582 


0.9253 


5.4 


1.39 


(1.02) 


SD170 


0.2733 


3 


8693 


0.9764 


5.0 


1.57 


(1.15) 


SH165 


0.2711 


3 


8578 


0.9260 


6.1 


1.36 


(1.03) 



poor abundance was preferred. Table [2] indicates that our 
metal-poor cases are no better than the other cases. The 
discrepancy might arise from the fact that their adopted ra- 
dius of HD 24712, log(J?/R0) x 0.231, is somewhat smaller 
than ours. 

The magnetic field strength for each best fit model is 
determined mainly by the requirement that the frequency 
difference between f8 and f6 be equal to the frequency dif- 
ference \i^{n4 + l;lm = 2) — u{n4, -\-l\lm = 3)1 which varies 
weakly as a function of Bp (Fig. [3| , where 714 is the radial 
order for the main frequency f4; i.e., 1^(714; Im ~ 2) =f4. (We 
note that at Bp = 0, \u{n;£ = 2) — v{n;l — 3)| is always 
much larger than the frequency separation between fS and 
f6 (or f7).) Since the magnetic field effect on pulsations is 
stronger in less dense atmospheres for a given frequency and 
Bp, the required magnetic field strength tends to be smaller 
for more massive or helium-depleted models. Considering 
the fact that the spectroscopically determined magnetic field 
strength is Bp « 4.4 kG (§[TI), a He-depleted model of AD165 
or AD170Z25 is better for the model of HD 24712. 

Fig. [5] shows the distribution of the mean deviations on 
the logi? — Bp plane for AD165 models, where R means 
stellar radius. Darker parts indicate smaller mean devia- 
tions. Theoretical frequencies are interpolated with respect 
to log R along the evolutionary model sequence, but not in- 
terpolated with respect to Bp (frequencies are obtained at 
every 0.1 kG). 

As stellar radius changes, the oscillation frequencies 
change with keeping Av approximately constant. Therefore, 
at a fixed Bp the mean deviation becomes small cyclically 
with varying log R, when regularly spaced observed frequen- 
cies (fml, fm2, and from fl to f6) become close to Im ~ I 
and 2 modes. This corresponds to the cyclic appearances of 
dark bands in Fig. [5] 

The dark bands are inclined because oscillation frequen- 
cies increase gradually as Bp increases. They are interrupted 
at Bp ^ 3.5 — 3kG, because oscillation frequencies jump 
there (see Fig. [S]). The darkness of the bands varies alterna- 
tively; i.e., the mean deviation is smaller when the frequen- 
cies fm2, f2, f4 and f6 are fitted with Im = 2 modes (as at 
Bp « 4.9 kG in Fig. [3} rather than with Im = 1 modes, be- 
cause in the former case the frequency f8 can be fitted well 
with an = 3 mode. 

Thus, we identify the main frequency f4 as an Im = 2 
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Shifted frequency (/xHz) COS 9 



Figure 6. Echelle diagram of frequencies for the best fit model in 
AD165 and ob served frequencies of HP 2471 2 from the WET cam- 
paign (pluses; iKurtz et aLll2005h and from iMkrtichian fc Hatzej 
1120051 ). Open symbols are model frequencies; circles for Im = 0, 
triangles for Im = 1, squares for = 2 and pentagons for Im = 3. 



Figure 7. The amplitude of radiative flux variation on the stel- 
lar surface as a function of cos for a model close to the best 
fit model of AD165, where 9 is the co-latitude measured from 
the magnetic axis. We have chosen one mode for each deg ree Im 
having a frequency close to one of the observed frequencies of 
HD 24712. 



mode because then the highest frequency f8 can be well fit- 
ted wit h an = 3 mode. T his identification is diff'erent from 
that of lKurtz et al.l (|2005l ): they identified f4 as a deformed 
dipole {Im = 1) mode mainly based on the consistency of the 
amplitude of rotational side-lobes. We will show below that 
due to a strong deformation of the amplitude over the stellar 
surface, the mean amplitude of the rotational side-lobes for 
f4 agrees with that expected from a, Im ~ 2 mode. 

Fig. |S] is an Echelle diagram for the model at 
(log i?, Bp) = (0.269, 4.9kG) in Fig. [5] It is the best from 
the family of AD 165 models for Bp > 3kG with the mean 
deviation of the model being 1.45 fiRz without Im ~ 
modes and 1.04 /iHz if Im =0 modes are included. (Al- 
though a smaller mean deviation of 1.21 /iHz is realized at 
{\ogR,Bp) = (0.261, 2. IkG), the magnetic field strength 
in the model is smaller than the observational estimates 
(§1), and the phase variation in the atmosphere is incon- 
sistent with observation.) Fig. [6] shows that although the 
large separation of the model agrees very well with obser- 
vations, the small spacing defined by = 1 and Im ~ "2 
modes is too large by about 3 fiHz. This discrepancy disap- 
pears if fl, 13 and f5 are fitted with Im = modes rather 
than Im ~ 1 modes. In this case, the nearly equal spacings, 
/5 - /4 /4 - /3 ^ /3 - /2 /2 - /I ^ 34/iHz, are 
realized by pairs of Zm = 2 and Im ~ 0. 

In previous investigations the identity of f8 was puzzling 
because it is separated from f7 (or f6) by only 17/iHz, half 
of the regular spacing. In our models /8 can be well fitted 
with an Im = 3 mode. It is not clear, however, why no other 
Im = 3 modes are detected. 



3.2 Amplitude modulation with rotation phase 

Observed pulsation amplitudes of roAp stars change with 
magnetic (rotation) phase in such a way that the ampli- 
tude is maximum at the phase of magnetic maximum. This 
is inte r prete d by the oblique pulsator model proposed by 
iKurtj jl982l ): i.e., an axisymmetric pulsation mode whose 
axis aligns with the magnetic axis which is in turn inclined to 
the rotation axis. As the star rotates one observes pulsations 
at varying aspect, which causes the amplitude modulation. 

HD 24712, like other roAp stars, shows an amplitude 
modulation with the amplitude maximum occurri ng at the 
magnetic maximum as found bv lKurtz et al.l (Il989l ) (see also 
.R yabchikova et al.llioOTi ') . Theoretically, amplitude modula- 
tion of a pulsation mode can be obtained by integrating the 
amplitude distribution over the stellar disc assuming vari- 
ous values of the angle between pulsation axis (i.e., magnetic 
axis) and the line-of-sight expected during a rotation period. 

Fig. [7] shows the amplitude of radiative fiux variation 
as a function of cos 6 for individual pulsation modes with 
various latitudinal degrees Im for the best model of AD165. 
6 is the co-latitude with respect to the magnetic axis. The 
frequency of each mode is close to one of the observed fre- 
quencies of HD 24712, though the amplitude distribution is 
insensitive to the pulsation frequency. Due to the presence of 
a strong magnetic field, the latitudinal distribution of pulsa- 
tion amplitude significantly deviates from any single Legen- 
dre function Pi{cos8). The amplitude around the equatorial 
region is strongly suppressed and it is more concentrated in 
the polar regions compared to the non-magnetic case. The 
Im ~ case is somewhat different from the other cases; it 
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Figure 8. Amplitude modulations with rotation phase in light 
variations predicted from the amplitude distributions on the stel- 
lar surface shown in F ig. [7] We have assum ed i = 137° and 
/3 = 150° according to iBagnulo et al.l lll995l '). and fi = 0.6 for 
the limb-darkening parameter. The quantity 7+ is defined as 
{A— I + A^i)/Ao with Ao and A±i being central amplitude and 
amplitudes of rotational side-lobes, respectively. 



has a broad peak around the equator, which significantly 
influences amplitude modulation as discussed below. 

It is remarkable that the amplitude distribution of the 
Ijn = 1 mode in the hemisphere is very close to that of the 
Im ~ 2 mode except that the latter (former) is symmet- 
ric (anti-symmetric) to the equatorial plane. Since we see 
mostly one magnetic hemisphere of HD 24712, we expect 
that amplitudes in light variations and rotational amplitude 
modulations of = 1 and = 2 modes are comparable to 
each other. 

Fig. |5] shows amplitude modulations predicted from the 
amplitude distributions shown in Fig. [T] assuming Bagnulo 
et al.'s (1995) parameters: i = 137°,/? = 150°. Magnetic 
maximum corresponds to the rotation phase of unity. All 
the modes except for Im = have maximum amplitudes at 
magnetic m aximum in agreement with the observations by 
iKurtz et al.l 0989). 

Despite a considerable difference in the amplitude dis- 
tributions between Im ~ 1 and Im = 3 (Fig. [7]), the ampli- 
tude modulation curves are very similar. This comes from 
the fact that for an odd mode contributions from compo- 
nents of ^ 55 3 to the light variation is very small com- 
pared with that of the 1 = 1 component, as discussed in 
ISaio fc Gautschvl l|2004l ). 

In contrast to the other cases, amplitude of the Im ~ 
mode is maximum at the magnetic minimum, although the 
modulation amplitude is small. This is caused by a broad 
peak of amplitude near the equator seen in Fig. [T] The prop- 
erty of the amplitude modulation of the Im ~ mode indi- 
cates that the largest amplitude mode of HD 24712 cannot 
be matched with a, Im ~ mode, although it is still possible 



that Im = modes are involved in the pulsations of this star 
without influencing the total signal very much. 

The observed ratio of the minimum to t he maximum 
ampli tude can be read as ~ 0.3 from Fig. 1 of iKurtz et al.l 
(Il989l ). which is consistent with the theoretical predictions 
for Zrri = 1, 2 and 3. More quantitative comparisons are pos- 
sible by using rotational side-lobes in the Fourier spectrum. 
The rotational side-lobes are characterized by the quantity 
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(4) 



where Ao refers to the amplitude of the central frequency 
and A-i and A+i refer to the amplitudes of the first lower- 
and higher-frequency side-lobes, respectively. (For our the- 
oretical side- lobes A^i — A+i because we do not include 
Coriolis force effects.) The value of 7+ for each mode in 
Fig. |8]is indicated in parentheses. The value of 7"*" depends 
on the latitudinal degree Im, but it is insensitive to the fre- 
quency; i.e., pulsation modes with the same Im but different 
radi al orders have simi lar values of 7'''. 

iKurtz et al.l (|2005| ) obtained 7"*" for each of the frequen- 
cies of HD 24712; among them the values for the three main 
frequencies f2, f3 and f4 are well determined with uncertain- 
ties less than 6%. The mean value from 2000 and 1986 data 
is 0.60 for f4, and 0.64 for f2. We have matched f4 and f2 
with Im ~ 2 modes which have 7"*" = 0.55 (Fig. |8]), not very 
different from the observed values. On the other hand, the 
mean value for f3 is 0.83, which is considerably higher than 
0.63 expected for an Zm = 1 mode. As discussed above and 
indicated in Fig. [S] f3 is also close to anlm ~ mode having 
7"'" = 0.99. The observed value of 7^ lies in the middle be- 
tween 7"'"s of = 1 and Im = modes. This could indicate 
the possibility for the frequency f3 to be a superposition of 
the Im ~ I and Z™, = modes. 

The amplitudes of the second rotational side-lobes are 
predicted to be about 10 times smaller than those of the 
first side-lobes. For example, the amplitude of the second 
side-lobes of f4 {Im = 2) is expected to be about 30 /xmag, 
which is smaller than the amplit ude of the highest noise 
peaks (80 /imag) in the analysis of iKurtz" et all l|2005h . 

Fig. [9] shows rotational modulations of pulsation phase 
at 5L — and at Kad = (changing from negative to posi- 
tive) for each mode shown in Fig. [T] where Vrad is the radial 
velocity (RV) at log r « —6. Except for the Im = case no 
rapid changes occur in pulsation phase as observed in some 
other stars (e.g. HR 3831 (HP 83368). [Kurtz fc Shibahashil 
ll986l : lBaldrv fc Beddin3l2000l ). This is because we see only 
one magnetic pole of HD 24712 with a configuration of 
(/3, i) = (150°, 137°). For the Im = mode rapid changes 
in pulsation phase occur at the rotation phases ~ 0.8 and 
~ 1.2. This comes from the presence of a nodal line at 
cos6' ~ 0.9 in the amplitude distribution (Fig. [7|. 

For the Im = 2 mode (the 2nd panel from the top in 
Fig. |5J, which corresponds to the main frequency f4, the 
phase of the luminosity variations is nearly constant. The 
phase of the RV variations is also practically constant in the 
rotational phase interval 0.75 — 1.25, although it changes 
gradually at other phases. Simultaneo us spectroscopic and 
photometric monitoring of HD 24712 (|Sachkov et al.ll2006l ) 
(see also Paper I) performed between rotational phases 0.87 
and 1.18 do not show phase variations, thus supporting theo- 
retical predictions. The phase of the luminosity variations is 
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Figure 9. Theoretical oscillation phase modulations in light and 
velocity variations for each mode shown in Fig. [7] In each panel 
oscillation phases at 5L = and at Vrad = (from negative to 
positive) are plotted as a function of the rotation phase, where 
Vjad refers to the radial velocity at log r ^ —6. 



Table 3. Phase lags in the outermost layers 



Line 



A(A) logrsooo 



Pr III 


5285 


-6.41 


0.44 


Pr III 


5300 


-6.35 


0.46 


Pr III 


6160 


-6.33 


0.42 


Pr III 


6090 


-6.30 


0.43 


Pr III 


6196 


-6.29 


0.43 


Nd III 


5294 


-6.03 


0.33 


Nd III 


5127 


-5.94 


0.31 


Nd III 


4927 


-5.79 


0.28 


Pr III 


6707 


-5.53 


0.37 


Pr III 


5844 


-5.51 


0.38 


Pr III 


5999 


-5.50 


0.38 


Nd III 


6145 


-5.49 


0.28 


Pr III 


6053 


-5.48 


0.39 


Nd III 


5410 


-5.44 


0.26 


Nd III 


5988 


-5.43 


0.27 


Nd III 


5845 


-5.32 


0.27 



always smaller than that of the RV variations, which means 
that the maxima (for example) of the luminosity variation 
precede the maxima of the RV variations. The maximum 
value of the phase difference is about 0.3 of the period and is 
realized when the rotation phase is around unity; i.e., when 
the pulsation amplitude is maximum (Fig. [8]). This means 
that velocity maxima lag luminosity maxima by about 30% 
of a pulsation period around the rotation phase when oscil- 
lati on amplitude is max i mum . 

iRvabchikova et all ()2007l 'l obtained phase lags of RV 
variations with respect to the hght variations for various 
spectral lines. These phase lags are always negative which 



means that luminosity maximum occu rs after the RV max - 
ima. For two ot her roAp stars. 10 Aql llSachkov et al.ll2008h 
and HD 101065 (|Mkrtichian et ahlbOO^ ). the observed phase 
lags are positive, in the range 0.4 — 0.6 (10 Aql) and 
0.16 — 0.19 (HD 101065), respectively. Allowing one cycle for 
the phase uncertainty, the negative phase lags in HD 24712 
may be converted into the positive ones. Table |3] lists some 
results for the lines formed in the most superficial layers, 
where optical depth of the formation of each line has bee n 
obtained by using a NLTE model of IShulvak et al.l l|2009t ). 
Generally the phase lag S(j) decreases with depth, indicat- 
ing oscillations propagates outward as discussed below. Ta- 
ble [3] shows that the phase lag around logr « —6 is about 
0.3 — 0.4, roughly consistent with the theoretical prediction 
for the main frequency. 



3.3 Amplitude and phase variations with depth 

IRvabchikova et al] (|2007l ) obtained amplitudes and phases 
of the RV variations in various spectral lines for the two high- 
est amplitude modes (f4 and f2) from time-resolved spectro- 
scopic observations of HD 24712. In the present paper we 
consider the results at a rotational phase of 0.94 (UVES 
observations close to the magnetic maximum). For a rep- 
resentative sample of Pril, Prill, Ndlland Ndlll lines, as 
well as for the Ha core, d epths of formation were calcul ated 
in NLTE approximation llMashonkina et al llioosl. |2009| ) us- 
ing a model atmosphere from lShulvak et al.l ( 20091 ) . In these 
iterative model calculations empirical stratifications of the 
chemical elements Si, Ca, Cr, Fe, Sr, Ba, Pr and Nd were de- 
rived for each iteration. NLTE calculations were performed 
for Pr and Nd only, because these elements are concentrated 
in the uppermost atmospheric layers, while other elements 
have a tendency to be accumulated close to photosphere 
where NLTE effects may be neglected. Final element distri- 
butions were then used for depth of formation calculations. 
The results have been converted to the relations between 
optical depth and RV amplitude/phase. For Yii lines no 
stratification analysis was performed, therefore, depth of line 
formation was calculated with a uniform yttrium distribu- 
tion using an yttrium abundance log(y/A'tot)=8.60 derived 
from abundance analysis. 

Figs l 101 aud i 111 compare RV amplitude/phase ~ logrsooo 
relations for f4 and f2, respectively, with the corresponding 
theoretical relations of Z™, = 2 modes for He-depleted best 
models (Table [5} . For Y ll and Fe ll lines with very low RV 
amplitudes pulsational phases are plotted twice: as defined 
originally in RyeijDchikQya et al. (2007), and shifted down- 
ward by a pulsation period taking into account the one- 
period uncertainty in the phase relation. Similar compar- 
isons for models without He-depletion are given in Figs ll2l 
and ll3l In comparing the observed and theoretical relations, 
the theoretical amplitude is multiplied by an arbitrary fac- 
tor, and the theoretical phase delay is fitted at the outer 
boundary (log r « —6) by applying an arbitrary shift. 

In the outermost layers oscillation phase gradually de- 
lays toward the outer boundary, indicating the presence 
of running waves propagating outward. This is consistent 
with the fact that all the observed oscillation frequencies of 
HD 24712 are larger than the acoustic cut-off frequencies of 
theoretical models. (If a reflective outer boundary condition 
is imposed, the phase in the outermost layers is nearly con- 
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Figure 10. Optical depth versus phase delay (top panel) and am- 
plitude (bottom panel) of RV variations for f4 (the main oscilla- 
tion mode) obtained from various spectral lines of HD 24712. For 
Y II and Fe II the phases shifted downward by a pulsation period 
are also shown. Various lines are theoretical relations for = 2 
modes having frequencies similar to f4 for some He-depleted mod- 
els (see Table [T] for parameters of each model). In calculating 
theoretical relations angles of (i,/3) = (137°, 150°) are assumed. 

slant.) The gradient of the phase variation in the outermost 
layers is weU reproduced by theoretical relations. We note 
that if the dS/dt term in equation ([2} is d ropped, the gra- 
dient becomes steeper, indicating that the lUnno fc Spiegell 
1 19661 ) theory is consistent with the observations. 

The phase distribution in HD 24712 has a jump between 
r K, 10"* and 10~^ for both f4 and f2, indicating the pres- 
ence of a quasi-node there. This property agrees with the 
theoretical phase variations for A-models with the T — r re- 
lation having a small temperature inversion at logr —3.5 
(Fig. [ij . A rapid phase change in S-models with the stan- 
dard T — T relation, however, occurs between r « 10~^ and 
10"'^, deeper than the phase jump derived from observa- 
tions. Thus, the position of the phase jump supports the 
T — T relation with a temperature in version at log r ~ —3.5 
as obtained bv lShulvak et al.l (|2009l ). 

Whether the phase of theoretical models jump down- 
ward or upward depends on the assumption of helium de- 
pletion, and sometimes on the strength of the magnetic field. 
Generally, phase tends to decrease inward at the quasi-node 
in helium-depleted models, and to increase inward in models 
without helium depletion. In the best model of the AH165 
series with Bp = 5.5 kG, the phase of RV variations for the 
main frequency f4 decreases steeply inward at logr ~ — 3, 
while it increases steeply at the same layer if Bp is slightly 
increased to 5.8 kG (Fig. ll2p . although for the secondary fre- 
quency f2 the phase increases steeply inward for both cases 
(Fig. I13p . However, such a strong dependence of the direc- 
tion of the phase jump on Bp does not occur in the other 
cases. 



Figure 11. The same as Fig. [TO]but for f2. 
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Figure 12. The same as Fig. llOl but compared with models with- 
out helium depletion. 



To obtain information on the helium depletion in the 
atmosphere of HD 24712, further observations are needed 
to fill the gap between log r — —4 and —3 . 

Towards deep interior, theoretical pulsation phase ap- 
proaches a constant; i.e., nearly a standing wave is realized 
in the deep atmosphere. However, the innermost Fe II data 
for the main pulsation f4 deviate considerably from the theo- 
retical curves. A possible cause may be related to the surface 
element distribution; Fe is concentrated around the equator, 
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Figure 13. The same as Fig.[T2]but for f2. 



while Y (as all rare earth elements , including Pr an d Nd) is 
concentrated near magnetic poles ()Liiftinger et al.l l2008). 

Oscillation amplitude increases rapidly in the outermost 
layers due to a rapid decrease in the gas density. In the lay- 
ers of —5 < log r < —4, theoretical amplitude variations 
roughly agree with the observed ones for both f4 and f2. In 
the outermost layers with log r < —5, however, theoreti- 
cal amplitudes deviate significantly from the observations. 
Observed amplitudes level off in the most superficial layers, 
the cause of which is not clear; it could be the effect of non- 
linear dissipation, or the density stratification could differ 
significantly from our simple models. 



4 CONCLUSIONS 

We discussed theoretical models for the oscillations of the 
roAp star HD 24712 (HR 1217). Observed frequencies are 
fitted well with theoretical ones for models whose positions 
on the HR diagram are close to the position of HD 24712 
determined by spectroscopy along with the Hipparcos paral- 
lax. The observed main frequency f4 is identified as a quasi- 
quadrupole {Im = 2) mode, whose amplitude on the sur- 
face is strongly concentrated in the polar regions and is 
suppressed significantly around the equator. The theoreti- 
cal amplitude distribution predicts a rotational modulation 
of the pulsation amplitude which is consistent with the ob- 
servations. 

We modelled for the first time the distributions of the 
phase and amplitude of RV variations as a function of atmo- 
spheric height and compared these with the observed distri- 
butions. The gradual outward increase of phase lag in the 
outermost layers is well reproduced by theoretical results ob- 
tained with a running-wave outer boundary condition. The 
presence of a steep phase change between log r « —4 and 



—3 favours a T — r relation with a temperature inversion at 
log ~ —3.5, rather than a standard T — r relation. 

Although our models agree with most of the observed 
properties of the oscillations in the atmosphere of HD 24712, 
we have recognized three problems to be solved in the future. 

(i) The observed amplitude of velocity variation levels off 
in the outermost layers with logr < —5, while in all the 
models we have calculated amplitude increases steeply out- 
ward in the layers with log r < —4. 

(ii) If the observed main frequencies are identified with 
Im ~ 2 and Im ~ 1 modes, to explain their nearly equal 
spacings the small frequency spacing defined as ^{n, Im = 
1) — 0.5[h'{n, Im ~ 2) + u(n —l,lm = 2)] must be as small as 
0.5 /iHz. For all the models calculated, however, the small 
spacing was always larger than ~ 3 /^Hz. The problem goes 
away if we assume the observed equal spacings of frequen- 
cies are produced by Zm = 2 and Zm = modes rather than 

= 2 and Im ~ 1 modes. However, this solution has a prob- 
lem; Im = modes have rotational amplitude modulations 
different from the modulation seen in the light variations of 
HD 24712. 

(iii) We have found no excited oscillation modes with 
frequencies appropriate for HD 24712; in other words 
all the modes examined are damped ones. High-order p- 
modes in roAp stars are generally thought to be ex- 
cited by the K-mechanism in the hy drogen ionization 
zone teal mforth et al.ll200ll : ICunhall200a ). However, the k- 
mechainsm does not seem to be strong enough to excite the 
super-critical high-order p-modes in HD 24712. It is inter- 
esting to note that the inefficiency of the K-mechanism is in 
common with the case of another cool roA p star HD 101065 
having many regularly spaced frequencies (jMkrtichian et all 
2008). Probably we need to find a new excitation mechanism 
for these coolest roAp stars. 
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